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Abstract 

We develop a new metric for quantifying end-to-end throughput in multihop wireless networks, 
which we term random access transport capacity, since the interference model presumes uncoordinated 
transmissions. The metric quantifies the average maximum rate of successful end-to-end transmissions, 
multiplied by the communication distance, and normalized by the network area. We show that a simple 
upper bound on this quantity is computable in closed-form in terms of key network parameters when 
the number of retransmissions is not restricted and the hops are assumed to be equally spaced on a 
line between the source and destination. We also derive the optimum number of hops and optimal 
per hop success probability and show that our result follows the well-known square root scaling law 
while providing exact expressions for the preconstants as well. Numerical results demonstrate that the 
upper bound is accurate for the purpose of determining the optimal hop count and success (or outage) 
probability. 

I. Introduction 

Determining the capacity of distributed wireless networks {i.e., ad hoc networks) is one of the 
most general and challenging open problems in information theory. Straightforward applications 
of known information theoretic tools and inequalities becomes intractable almost immediately 
and have hence yielded little in the way of results. This motivates the exploration of approaches 
to describing ad hoc network throughput that, while falling short of strict information theory 
upper bound standards, do provide insight into the fundamental trends on achievable throughput. 
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appeared in part at the 2009 Allerton conference. Manuscript date: September 30, 2009. 
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A. Motivation and Related Work 

The main line of present inquiry is to consider the transport capacity of an ad hoc network, 
which quantifies the bits per second that can be reliably communicated over some distance in 
the network [1]. A significant merit of this approach is that it describes important aspects of the 
capacity region of all possible rate pairs in the network - notably how the region scales with 
the number of nodes n - and has provided high-level insight on how different network scenarios 
and approaches may affect the scaling law. Notwithstanding special cases that suggest more 
optimistic scalings [2], [3], [4], it is generally agreed that the scaling law in ad hoc networks 
is 0(1/71) and can be achieved with nearest-neighbor routing through the network [5], [6], [7], 
[8], [9]. Despite this considerable progress, a common limitation of this line of research is that 
the main results are typically asymptotic scaling laws that do not easily allow capacity tradeoffs 
to be made based on the salient network parameters and design choices. Although many efforts 
have been made to better understand the transport capacity preconstants, e.g. [10], [11], [12], in 
general tractability has suffered and/or several essential aspects of ad hoc networks have had to 
be abstracted out. 

An alternative approach pioneered by the current authors and others has focused on computing 
achievable rate regions and network densities for different architectures and technologies. A key 
to this approach is to assume that the interferer locations are random, in particular that they are 
Poisson distributed over the plane. Although this is an idealization, it accurately models the sce- 
nario where transmitters are randomly scattered and uncoordinated, which seems to be quite a bit 
more realistic than other popular alternatives, such as assuming they are placed deterministically 
on a regular grid. This approach has allowed closed-form derivation of outage probability as 
well as the maximum allowable transmission density at a specified outage probability and data 
rate, the latter of which we have termed the transmission capacity [13], [14], [15]. The analytical 
tractability of this approach (using tools from stochastic geometry [16], [17]) has to date allowed 
quantitative design tradeoffs for spread spectrum [13], [18], [19], interference cancellation [20], 
[21], multiple- antenna architectures [22], [23], [24], [25], cognitive radio and capacity overlays 
[26], [27], [28] and power control [14], [29]. 

A common shortcoming of the transmission capacity approach is that it considers a typical 
network snapshot, and hence only simultaneous single hop transmissions. The primary goal of 
this paper is to address this limitation, and we develop a related metric which we term the 
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random access transport capacity because it presumes packets are transported end-to-end over 
some distance R, but assumes independent locations and transmissions for the interferers. In 
that sense, one contribution of the paper is to find a middle ground between the transport and 
transmission capacity approaches into the capacity of multihop wireless networks. As one might 
expect, the results of this paper follow the Q(y/n) scaling law, but provide additional information 
on the preconstants. 

Complementary to the transport capacity research, some recent work has formulated the 
multihop capacity problem as a line network without additional network interference [30], [31]. 
Both of these papers agree that numerous hops are helpful only in the "power-limited" regime, 
that is where the spectral efficiency is low and overcoming noise is the primary concern. Both 
also find that in the "bandwidth-limited regime" - when the SNR is high - that additional hops 
decrease the end-to-end throughput due to the use of extra time-slots. Two notable theoretical 
results are i) that the end-to-end capacity scales as O(logM) for M hops [30] and ii) that by 
using optimal time-sharing on each hop m — 1, . . . , M to achieve the per-hop capacity c m , the 
end-to-end capacity with no interference increases from C no int = min m (c m ) in [30] to 



as derived in [31]. Intuitively, ([T]) tells us that each additional hop lowers the capacity unless 
the hop capacities c m increase by a sufficient amount (due to higher SNR). Hence this equation 
also captures a basic tradeoff between higher per hop capacity (more short hops) and the desire 
for fewer total dimensions or transmission (fewer, and hence longer, hops). This tradeoff has 
been considered from many possible perspectives in [32], but in the present paper we are able 
to quantify it exactly in terms of a transport capacity-like metric, and with interference. 

Finally, we note that a recent paper [33] has considered multihop capacity in a Poisson field 
of interference. The key distinction in this work is that we focus on achievable end-to-end 
throughput and optimal transmission strategies and hop count rather than the stability of queues. 
As [33] argues, the delay-optimizing number of hops appears to also be throughput-optimizing, 
so our conclusions are not contradictory. We also consider noise whereas they do not, which is 
quite important in the power-limited regime where multihop is beneficial. 




(1) 
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B. Contributions and Organization 

In this paper we develop a new, quite general model for end-to-end throughput in a multihop 
wireless network. We term the resulting metric random access transport capacity since the 
analysis requires all transmissions to be independent, which precludes cooperative transmission 
scheduling among the nodes, since this would generally couple transmissions and the active 
transmitter locations would no longer be independent. Note, however, that the model does not 
preclude cooperative or multipacket reception, although we do not consider such approaches in 
this paper. The general model includes arbitrary paths of M hops and an end-to-end delay/energy 
constraint quantified as the total allowable number of transmissions per packet, A. A closed-form 
upper bound on random access transport capacity is obtained by making the M hops equidistant 
and by letting A — > oo. This upper bound is tight for moderate values of A. 

From the upper bound, the optimal number of hops M* can be computed as a function of the 
network parameters - in fact this was a prerequisite to finding the upper bound. For example, 
we can show that M* oc Rf3^\f\, where R is the end-to-end distance, f3 is the per-hop SINR 
required for successful communication at the transmitted rate, a is the path loss exponent, and 
A is the spatial density of transmitters. Additionally, the optimal per-hop success probability (or 
equivalently, outage probability) and optimal average number of retransmissions per hop can be 
directly found. This exposes an optimal tradeoff between highly reliable transmissions - achieved 
with short hops (large M) and hence a low per-hop outage - and unreliable transmissions, which 
make up for the retransmissions by requiring fewer hops. The upper bound is shown to be accurate 
for the purposes of predicting the optimal number of hops M* and the optimal per-hop success 
(or outage) probability. 

An advantage of the novel framework developed in this paper is that the end-to-end throughput 
(here, the random access transport capacity) uses a small, finite, and computable number of 
hops, and results are given as a function of the S — V separation R. We prove that the random 
access transport capacity follows the square-root scaling law, while providing easily computable 
preconstants. We expect that this approach can be extended to better understand and compare 
many specific candidate technologies and design approaches in multihop wireless networks. 

The rest of this paper is organized as follows. First in Section [TTJ, we define the overarching 
network model and metrics used in this paper. We then develop and derive the upper bound in 



Section III which is the main result of the paper. Implications and observations based on the 
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upper bound derivation are discussed in Section IV Numerical results are given in Section |V} 



in which the upper bound is evaluated against a numerical exact solution (for finite A) and the 
above insights from the results are discussed. 

II. Network Model and Metrics 

We consider a large wireless ad hoc network where a typical source node S drawn from a 
homogenous 2-D Poisson point process (PPP) of intensity A wishes to transmit to a destination 
node V that is a distance R away in a random direction and is not part of the PPP All transmitters 
(desired and interfering) have fixed transmit power p, which in order to simplify notation is more 
precisely the average radiated power at a distance of 1 meter from the transmitted The noise 
power is r] and the channel strength is determined by path loss and fading, so the received 
power at a distance R is p\R~ a , where a > 2 is the path loss exponent and x results from iid 
Rayleigh fading, which means \ ~ exp(l). The SNR is defined in this paper as SNR := pR~ a ji] 
regardless of the number of hops taken, so dividing the route into hops improves the per-hop 
SNR, as it should. 

A. Single Hop, Single Transmission 

The simplest case is to consider a single transmission over the entire S — V distance R. In 
this case, the received SINR, where the reference receiver is defined to be at the origin, is 

SINR = - PXoR v " . (2) 

The interfering node locations in (|2]) are drawn from a stationary (Poisson point process (PPP) 
on the plane of intensity A, denoted II = {Xj}, where each X 4 E M 2 is the location of an 
interfering transmitter. In the single-hop model, all interferers are sources, themselves. This is a 
realistic model assuming the transmitting nodes in the network are randomly and independently 
located and do not cooperate. 

Under this model, the probability of success for sending a packet from source to destination 
can be found to be 

p s = p s (l) = exp - A/^ Q i? 2 } , (3) 

'The authors have proven in [14] and [29] that pairwise power control does not have a significant effect on outage probability 
or transmission capacity. Path loss models that avoid the singularity for R — > are discussed in [34]. 
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where K a = 2n 2 /(a sin(27r/a)). Obtaining this expression is crucial to the results in this paper 
but nontrivial: for the derivation the reader is referred to [35]. Note that a similar exact result 
exists for Nakagami fading [22] as well as a tight bound on success probability with path loss 
only (no fading) [13], but we will use the expression for Rayleigh fading in (|3]) throughout this 
paper. 

This single hop model has been fairly well-studied in recent years. For such a model, we 
define the random access transport capacity, in this case for single hop, to be 

C sh = Ps Alog(l + P)R. (4) 

In words, C sh gives the density of successful transmissions at rate log(l+/3) that can span a S—V 
distance of R: therefore C has units of bps/Hz/m. As R increases, p s — > exponentially fast per 
(|3]). If p s is fixed to be 1 — e for some maximum outage probability e, the R term is dropped in 
Q, and ^ is inverted to find A, then Q gives the so-called transmission capacity. The goal in 
this paper is to move beyond the single-hop, single-transmission model to a network that allows 
multiple hops and multiple transmissions per hop, while retaining some of the tractability of the 
transmission capacity model. 

B. Multiple Hops, Multiple Transmissions per Hop 

Now, allow the S — V distance R to be subdivided into M hops having not necessarily 
equal distances r m , where R = J2m=i r m- On each hop, retransmissions are allowed and 
acknowledgments are used. Therefore, the number of transmissions on each hop m is denoted 
by the geometric random variables (Ti(M), . . . , T M (M)), where each T m (M) E {1, 2, . . .}, and 
the total number of transmissions required to move a packet from S to V is 

M 

T{M) = Y,T m {M). (5) 

m=l 

The notion of outage comes in naturally by constraining the total number of transmissions per 
packet, which is functionally a constraint on the total amount of delay and energy per packet, for 
example due to a timeout. If the total transmission per packet allowance is A then T(M) < A is 
required for successful transmission, so an outage event occurs when T(M) > A. Therefore, the 
total number of actual transmissions per packet is min(T(M), A), or more compactly, T(M)AA. 
We assume that each S—V pair only transmits a single packet at a time along the entire multihop 
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path: in other words there is no intra-route spatial reuseQ The interference density is still A as 
in the single-hop case, but each packet is sent T(M) A A times, causing the effective rate per 
S — V pair to be log(l + j3)/(T{M) A A) to account for the retransmissions. In other words, 
if the transmitted bit rate at S is log(l + (3), because S has to wait T(M) A A packets before 
transmitting again, the effective S — V data rate is lowered by this factor. This also presumes 
that S is instantaneously notified when the packet is received successfully at V so it can send 
the next packet. 

Definition 1: Random access transport capacity. The random access transport capacity C 
for a multihop wireless network is the maximum average source to destination rate that can be 
sustained reliably over a distance R with at most A transmission attempts per packet, normalized 
by the area of the network. Formally, it can be characterized by maximizing over the number 
of hops as: 

C(A) = max P(T(M) < A) ■ X ■ JS^-P^t R (6) 
v ' m 6 {i,..,4} ~ E[T(M) A A] 

= Alog(l + /3)i? max = A } (bps/Hz/m). (7) 

6V HJ Me{i,...,A] E[T(M) A A] w 

This quantity therefore determines the end-to-end rate that can be supported when the source 
nodes have density A. The challenge is in computing the probability of success P(T(M) < A) 
and in computing the expectation of a minimum, i.e., E[T(M) A A). In the general model thus es- 
poused, T is the sum of M independent but not identical random variables (Ti(M), . . . , Tm{M)), 
where T m (M) ~ Geo(p s (M)), for p s (M) to be defined in ([9]). For each transmission per hop to 
have independent success probability requires sufficient diversity in the interference and signal 
strength per transmission attempt, which could possibly be achieved through diversity techniques 
such as frequency hopping. Even with these simplifications, the general model does not appear 
to be tractable. In the next section, we develop an upper bound on the random access transport 

2 This is an important simplification, as it avoids dependent simultaneous transmissions along a path, as well as difficulties 
in accounting for non-Poisson interference. From a network-wide perspective the lack of spatial reuse for a given path does 
not affect the transport capacity because other S — T> pairs can utilize the space instead. Strictly speaking, the relays must be 
deterministically placed on the line between S and T> rather than drawn from If in order to maintain a PPP for the interference. 
Since exactly one node transmits per route at a time, the interferers in each time slot form a PPP in this case. We also note that 
the spatio-temporal correlations in a Poisson field of interference have recently been studied in [36]: here, the effect would be 
to increase the independence as M increases. 
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capacity that is computable in closed-form. 

III. Upper Bound: Guaranteed End-to-End Delivery, No Delay-Energy 

Constraint 

A. Establishing an upper bound 

To develop an upper bound on the random access transport capacity, we make two important 
simplifications. First, we assume the hop lengths are all equidistant, that is r m = R/M, for 
m = 1, . . . , M. This is best-case because p s in (|3]) is o(e~' r ™), so the increase in p s on short 
hops is outweighed by the decrease in p s on longer hops. In other words, any perturbation from 
equal-length hops increases the end-to-end outage probability, all else being equal. 

This assumption allows much improved tractability because now the probability of success 
is the same for all hops, so the number of transmissions T m (M) required on a given hop 
for success becomes iid geometric with parameter p s (M) to be defined in (|9]). The second 
important simplification is that we relax the delay-energy constraint. Formally, we let A — > oo in 
P(T(M) < A) /E[T(M)AA]. This simplifies both the numerator and denominator: the probability 
of end-to-end success P(T(M) < A) tends to 1, and E[T(M) A A] simplifies to E[T(M)]. These 
considerations yield the following Lemma and Corollary. 

Lemma 1: For all A e Z+ and all Me {1, . . . , A}: 

P(T(M) < A) 1 

E[T(M) A A] ~ E[T(M)]' ( ) 

Proof: See Appendix. ■ 

Although the numerator P(T(M) < A) -> 1 and the denominator E[T(M) A A] -> E[T(M)] 

both increase, the numerator does so slightly more quickly, resulting in an upper bound. Under 

the upper bound assumption of equally spaced hops (r m = R/M), the per-hop probability of 

success becomes: 

p s {M) := exp |-^M" Q - \R 2 f3«K a M~ 2 ^ , (9) 

and the probability mass function of the number of transmission attempts T m (M) on each hop 
m is geometric with probability of success p s (M): 

P(T m (M) =t) = (l- p s (M))*-V s (M). (10) 

The average number of attempts per hop is therefore E[T m (M)] = l/p s (M). Because the M hops 
are iid it follows that the expected number of total transmissions required to move a packet from 
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source to destination is E[T(M)] = ME[T m (M)] = M/p s (M). Combining these observations 
yields the following upper bound on the capacity. 
Corollary 1: An upper bound on the capacity (|7]) is 

C ub (A) = XlogCl + /3)R max (H) 

^ J &V HJ MG{1,..,A} M 

Note that although we relax the A constraint in Lemma [T] inside the optimization in |7]) to 
establish the upper bound, we retain a finite A in the range of M G {1, . . . , A}. 

5. Optimal number of hops 

The next step is to determine asymptotic optimal number of hops M* that maximizes C nb (A) 
as A — > oo. This is given by the following theorem and corollaries, which give a positive 
real value for M* although the actual quantity would necessarily be an integer. We leave it 
as a continuous quantity in this paper for tractability and generality. Once M* is computed, 
the optimal average number of transmissions required per hop can be readily determined as 
E[T m (M*)] =p~\M*). 

Theorem 1: Optimal number of hops, M*. The number of hops M* that optimizes the asymp- 
totic transport capacity density upper bound lim^oo C ub (A), i.e., 

M* := arg max ^AMl (12) 

is the solution to the equation 

M a - 2\f3 2/a K a R 2 M a ~ 2 - ^-a = 0. (13) 

P 

This results in closed-form expressions^] for M only when a G {3, 4, 6, 8}, in which case M* is 



the largest positive root of ([T3J. 
Proof: The objective is 

= J- exp(-A; 1 M- Q - k 2 M~ 2 ) (14) 

where k\ = = and k 2 = X(3 2 ^ a K a R 2 . Taking the derivative and setting equal to zero 
eventually gives 

exp(A;iM~ a + k 2 M' 2 ) (l - k x aM~ a - 2k 2 M~ 2 ) = (15) 



3 By closed-form, we mean direct computation is possible using only basic arithmetic operations, simple trigonometric 
functions, and radicals. 
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which gives that 



1 - haM- a - 2k 2 M~ 2 = 



(16) 



or in polynomial form 



AT - 2k 2 M a ~ 2 - ha = 0. 



(17) 



By the Abel-Ruffini theorem, a formula solution to a polynomial equation only exists for when 
the degree of the polynomial is 4 or less. The path loss exponent a is of physical interest only 
when 1.5 ;$ a ^ 6. In this range, M can be found in closed-form only for a E {2,3,4,6}, 
although it can also be found in principle for a E {1,8}. The solutions for a = 6 and a = 8 
would follow similarly to the a = 3 and a = 4 solutions. 



The corollaries for a = 3 and a = 4 follow immediately. A solution in terms of K a can 
be found for a = 2, but for a — > 2, K a — » oo so this result does not exist. Intuitively, when 
a < 2 the expected interference at any point in the network is infinite since the expected number 
of interferers increases as d 2 but the interference fall off is d a . Thus, the optimal number of 
hops for a = 2 is technically infinite. In the sequel we write M* = M*(a) to emphasize the 
sensitivity of the optimal number of hops on the path loss exponent. 



Corollary 2: Solving ( 13 ) with a = 3 yields two possible solutions, depending on the polarity 



of the equation's discriminant D — -± — ^f. If D > 0, 



where 



AT (3) = p*R 



f 



2p 



2p 




3r]\ 2 8Ki 



27 



A 3 , 



and K 3 = w 7.6. When D < 



AT (3) = 2 V2Aif 3 ^|^ cog (1 os 
3s \ 3 



(18) 



(19) 



3\/3r/ 



4V2p{XK 3 



(20) 



Proof: See Appendix. 



4 That is, there exist empirical studies that have found a to be in this range 
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Although the expressions ( [18] ) and ( j20] ) at first appear to be quite different, in fact they have 
a quite similar dependence in terms of the main parameters of interest. Finally, we have the 
following corollary for a = 4. 



Corollary 3: Solving ( |T3| ) with a = 4 yields a single maximum positive, real solution for M* 
which is 



M*(4) = p<R 



\ 



7T Z 



A Y + \/ A2 4 



7T 4 477 



(21) 



since K 4 = n 2 /2. 



Proof: See Appendix. 



C. The upper bound, C nh 

The hop-optimized random access transport capacity can now be given in closed-form for any 
choice of system parameters since 



mb Alog(l + /3)i? _ X1 , mD p s (M* 



C ub = F ° v — = Alog(l + /5) J R^ (22) 

and we now have found exact expressions for p s (M*) and M*. Although a few assumptions are 

made to get to closed form - fixed equidistant relays, randomly located interferers, independent 

retransmissions, all transmissions at the same rate - this allows a simple closed-form end-to-end 

expression for multihop network throughput with both noise and interference. 

For example, consider a = 4, in which case at high SNR M* — > fltRit^/X, which yields with 

(|9]) the following quite simple exact expression for transport capacity 

VAlog(l+/5) f V 1\ 

ex P -— JT^-o • ( 23 ) 



lim C uh * 

SNR-*oo 



a=4 



7T/3S "V P^ 4 A 2 2 



IV. Upper Bound Implications 

A. Initial Insights from Theorem 1 

It can be immediately observed that all the derived solutions for M* are precisely proportional 
to R/3*. This makes sense: the optimal number of hops should scale linearly with the source- 
destination distance since this keeps the per-hop distance and hence the per-hop SNR and SINR 
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constant for a certain p, 77, A. Similarly, since the received SINR is proportional to the received 
(desired) power P r over a distance d, we can see that 

fM\ a 

P r=P d =P[-r) ■ OA) 

Since P r must scale linearly with (5 with the interference (A) held constant, 

f M\ a 

p[ — J = constant • (3 M oc /?«. (25) 



Less obviously, there is a trend in the expressions that M oc \/X, which also makes intuitive 
sense. As the density increases, the interferers' relative distance (statistically) decreases as V% 
requiring a shorter communication distance by the same amount in order to maintain the same 
SINR. Hence, R/M oc A~ 2 or M oc y/X. This trend is more difficult to see in the a = 3 case 
for D < ([20]) but since 

cos ( -arc cos x I ~ ^—^ + I 1 — - I x, (26) 




3 J 2 

one of the two terms has V% the other having a A -1 term. As A increases, the A -1 would 
decrease leaving M oc VX 

In summary, from Theorem [T] and its corollaries, one can immediately determine the optimum 
number of hops for any plausible integer path loss exponent in terms of the salient network pa- 
rameters. Non-integer values could be accurately estimated by standard interpolation techniques. 

B. Optimal success probability 

Armed with the optimal number of hops M* , it is straightforward to determine the optimum 
success probability per hop, p s (M*). Interestingly, it can be expressed in two ways, where the 
first expression depends on the A but not SNR (neither noise nor interference power affect this 
expression), and for the second expression, vice versa. 

Lemma 2: The optimal success probability per hop p s (M*) as a function of M* is 

Ps (M*) = (27) 
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Proof: Since the optimal number of hops M* should satisfy ( [17] ), we can write the following 
two expressions 

1 °'c 2 i 



kxM' 



1 - 2koM~ 2 



a 

1 - fcjM" 



(29) 
(30) 



Plugging (29) into @ gives (27), while plugging (30) into ph gives (28). 



The two expressions for p s (M*) can be seen as the p s which optimally balances between inter- 



ference ( [27] ) and noise ([28]). Note additionally that the optimal number of average retransmissions 
per hop is given by T m (M*) = l/p s (M*), which is the same for each hop. 

C. Relation to general transport capacity: shared scaling law 

Although the primary merit of the approach of this paper is the ability to compute end-to-end 
rate results in closed-form, the results here also obey the well-known Q(^/n) scaling law for 
transport capacity [1], [8], or in the lexicon of this paper, 0(vT). This fact is formalized in the 
following theorem. 

Theorem 2: The random access transport capacity upper bound is 0(\/A) and can be stated 
formally as 



lim 



C ub e-5log(l+/3) 



(3D 



gives 



Proof: First, it can be shown that M* = c ^/\ for some constant c . Rearranging ( flT] ) 

(3a 



2\(3«K a R z M- 



(32) 



SNRM" 

For this equation to hold as A — » oo, clearly M — » oo, which means that the quantity on the left 
of ( [32] ) must approach 0. Therefore, M = a/A • \j2K a f3^R, i.e., c = y / 2K a (3^R. 

Second, from ( [28] ) it can be observed that p s — > exp(— |) as A — >• oo, whereas ( [27] ) can be 



used to verify the value for c . Using (22) and the scaling results just derived gives 



lim C ub 

A^oo 



lim 

A— >oo 



\p s (M*)log(l + P)R 
M* 



Aexp(-|)log(l + /3) J R 



v / Ae-hog(l + /?) 



(33) 
(34) 

(35) 
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This theorem suggests that the random access transport capacity even subject to the assump- 
tions we have adopted in this paper is order optimal: it follows the square-root scaling law for 
source-destination transmissions in large wireless networks. That is, unscheduled, channel-blind 
transmissions can achieve - given well-positioned relays - a transport capacity that scales the 
same as optimally scheduled nearest neighbor routing. This complements a similar conclusion 
for slotted ALOHA in regular networks [12]. 

V. Numerical Results and Insights 

We set R = 1, (5 = 3, and as a default a = 3, while the other parameters are varied. Note 
that the units of R and A are arbitrary, but in a circle with S and V on opposite sides separated 
by a distance R, there would be on average \-k(R/2) 2 interferers. Because only together do 
the values of A and R provide a description of the interference environment, A ^> 1 may be 
reasonable to model a dense network or a long S — V path when R — 1. 

Fig. [T] demonstrates that the random access transport capacity upper bound is tight for a small 
and large number of allowed attempts A, but fairly loose in between. This is perhaps to be 
expected since the bound presumes A — > oo. Figs [2] and [3] show the upper bound vs. the actual 
random access transport capacity (computed numerically) for A = 6 and A = 12, respectively, 
vs. the number of hops used, M. As would be expected by Fig. [T} the former is loose while the 
latter is fairly tight. What is interesting is that in both cases the maximizing M is very similar. 
In other words, even when the capacity upper bound is loose, the value of M* that it predicts 
appears to be quite accurate. This is confirmed in Fig. |4j the M* from the UB and the actual M* 
are always within one hop of each other and asymptotically the same. Interestingly, the actual 
M* is not necessarily monotonic in A, as Fig. ^demonstrates. 

Now focusing on the upper bound, first consider the result of Theorem 1, which gave the 
optimum number of hops M*. In Fig. [5] we observe that the M* is indeed 0(A/A) and that it 
obtains this scaling in the A G (.1,1) range over a wide range of SNR and for both a = 3 
and a = 4. The scaling law here kicks in faster at high SNR, because at low SNR a higher 
contention density is required before the network becomes interference limited. Although the 
path loss exponent does not appear to affect the convergence, smaller path loss exponents require 
more hops to be taken in order to mitigate the cumulative interference. This may seem counter- 
intuitive because larger path loss exponents cause faster attenuation of the desired signal, which 
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might seem to imply that more hops would be needed to maintain the signal strength, but this 
is not the case, even at very small interference densities. 

In Fig. [6] we plot the random access transport capacity upper bound where we do not optimize 
for M , but rather see how it varies with the number of hops. The optimum throughput occurs 
at the integer value of M closest to M*, as expected. The random access transport capacity is 
considerably higher overall with larger A (i.e., source density) and SNR, as would be expected, 
and achieving that higher end-to-end throughput requires significantly more hops in order to 
maintain a high enough SINR on each hop. On the contrary, if a small number of hops was 
used in this case, despite the high SNR, the end-to-end throughput would be much lower than 
in the lightly loaded (small A) case. Optimizing the number of hops - or at least getting close 
to M* - is very important. 

Turning our attention to Theorem 2, we observe in Fig. [7] that the scaling law of Q(y/X) holds 
for large A. Similar to Fig. [5] the scaling law kicks in faster at high SNR and the throughput 
increases with a. Put optimistically, this means at low SNR (or low interference density) that the 
throughput actually increases quite a bit more quickly with A than the scaling law would suggest. 
Until the network is saturated with interference, more users can be added without having their 
throughput decrease. 

Appendix 
Proof of Lemma 1 

Fix a positive integer A and p 6 [0, 1]. Define T(M) = T M ~ Pascal(M, p) as the number 
of trials required for M successes, and S A ~ Bin(A,p) as the number of successes in A trials. 
To prove the lemma it suffices to show f(M) > for all M — 0, . . . , A, where 

f(M) := pE[T M A A]- MF(T M < A). (36) 

The inequality is trivially true for M > A and M — 0, and is straightforward to show for M = 1 
and M = A. We first prove some key relationships for binomial random variables. 

A 

P Y j nSn = M) = F(S A+1 >M + 1), M = 0,...,A-1 (37) 

n=M 

F(S A <M-1) = (1 - p)F(S A -i < M - 1) + pF(S A -i < M - 2) (38) 
(l-p)(M + 1)F(S A = M + 1) = p(A — M)F(S A = M) (39) 
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To see ([37]), use the fact that pP(S n = M) = ¥(T M+1 = n + 1) and {T M < A} = {S A > M} 
to show: 

A A 

p p (S«. = M)=J2 P ( T m+i = n + 1) = P(T M+1 <A + 1)= F(S A+l > M + 1). (40) 



n=M 



n=M 



To see ( |38| ) simply observe the equality of the events, where X A indicates success on hop A: 
{S A <M-1} = {(S A -i < M - 1) n (X A = 0) U (S A -i < M - 2) n (X A = 1)}. (41) 



To see (f39j) use the "mixed product" identity: (M + 1) ( M ^J = (A — M) (jj) . That is, with D 
as the difference of the two sides in ( f39| ): 

£> = (l-p)(M + l)F(S A = M + l)-p(A-M)F(S A = M) 

= (1 -P)(M + 1)( M ^ jV+^l -p)^ 1 -p(A-M)(£V(l -p) A " M 



P 



M+l. 



1-p) 



A— M 



(M + l) 



M + l 



- (A - M) 



M 







(42) 



Using ( [37] ) we have the following development for f(M): 
f(M) = pE[T M AA]-MP(T M < A) 

/; 1 ^ nP(T A/ = n) + A¥(T M > A) ) - MP(T M < A) 



Vn=M 
.1 



nF(S n -x = M - 1) + pA(l - P(T M < A)) - MP(T M < A) 



n=M 



pA + M P J2 F (S n = M) - (pA + M)P(T M < A) 

n=M 

pA + MF{S A +i > M + 1) — (pA + M)P(S A > M) 
(p4 + M)¥(S A < M - 1) - MP(S A+1 < M) 



(43) 



We will prove /(M) > by showing that A(M) = /(M + 1) - f(M) > for all M = 
0, . . . , A - 1. Combined with /(0) = 0, this immediately yields that f(M) > for all M = 
0, . . . , A, establishing the upper bound. Consider the difference A(M) ee /(M + 1) — f(M) and 
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apply ( |38] ) twice: 

A(M) = [{pA + M + 1)¥(S A < M) - (M + l)¥(3 A +i < M + 1)] - 
[{pA + M)F(S A < M — 1) - MF(S A+1 < M)] 
= [(pA + M + l)¥(S A < M) - (M + 1) ((1 - p)F(S A < M + 1) + P F(S A < M))] - 

[(pA + M)F(S A < M - 1) - M ((1 - p)F(S A < M) + P F(S A < M — 1))] 
= ((1 - p)(2M + 1) + P A)F(S A < M) - ((1 -p)M + P A)¥(S A < M - 1) - 

((l-p)(M + l))P(^<M + l) (44) 
Now use the easy facts that: 

F(S A < M - 1) = P(S4 < M) — P(S A = M) 

P(5 A <M + 1) = F(S A < M) +F(S A = M + 1). (45) 
Substitution of ( |45] ) gives: 

A(M) = ((1 -p)(2M+ 1) +pA)F(S A < M) - 

((1 - p)M + pA) \P(S A <M)- F(S a = M)] - 
((1 -p)(M+ 1)) pP(S A < M) + P(5 A = M + 1)] 
= ((1 — p)M + pA)F(S A = M) — (1 — p)(M + 1)F(S A — M +1) (46) 
Apply (f39|) to the previous equation 



A(M) = ((1 — p)M + pA)F(S A = M) — (1 — p){M + 1)F(S A = M + 1) 

= ((1 - p)M + pA)P(«SU = M) - (1 - p)(M + V ^^mI^ a = M ) 

= F(S A = M)(((l-p)M + pA) -p(A- M)) 

= ¥(S A = M)(M -pM + pA-pA + pM) 

= MF(S A = M) > (47) 

Proof of Corollary 2 

For the D > case, a necessary condition on A and the transmit power can be formulated by 
inserting k x = -J^ = g|g and k 2 = X(3 2/a K a R 2 to get 
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The roots of such a 3rd order polynomial consist of 1 real and 2 complex roots. The real root 
is given by 

(49) 



which results in (18). 



When D < there exist three distinct real roots. The necessary conditions (48) are simply the 



converse of ( |48| ), but the polynomial solution requires some trigonometric calculations. Define 

3Jfei 



2 /2Jfe2 



2C 



arc cos 



3*4 

2C 



then the three possible solutions are 

yx = 2(3 cos 



y-2 



2/3 



2(3 COS ( + yj 



2(s cos 



47T 



(50) 

(51) 
(52) 
(53) 



For these three real roots, Vieta's Theorem states that for a 3rd order polynomial a 3 y 3 + a 2 y 2 + 
a\y + a that the roots will satisfy y\yiy?, = — —, which gives from (TT7b that 



V1V2V3 = 3A:i ^> 2/i2/2?/3 > 0. 



(54) 



Thus the roots must either be all positive or comprise one positive and two negative. They cannot 
all be positive since 2/12/2 + 2/i2/3 + 2/22/3 = 0' a l so a consequence of Vieta's Theorem. Hence, there 
is a single positive root, which is yi, and by expanding ( and <p can be expressed as M* = yi, 



given in (20). 



Proof of Corollary 3 
The polynomial for a = 4 becomes 



M 4 - 2k 2 M 2 - 4fci 



0. 



(55) 



Using a dummy variable for M 2 , this can be easily solved using the quadratic formula to give 



M*(4) = ±p*R 



\ 



4r/ 



(56) 



of which only the solution given in ( |2Tj ) is both positive and real 
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TABLE I 

Summary of notation and parameters 



A 


interference density, intensity of PPP II 


a 


path loss exponent (a > 2) 


K a 


A constant, K a = 2-k 2 /{a sin(27r/a)) 


13 


target (required) SINR per hop 


R 


S — T> transmit distance 




transmit distance on hop m 


p 


transmit power 


n 


noise power 


SNR 


S — V end-to-end signal to noise ratio, SNR := pR~ a /ri 


Ps(M) 


probability of success over a single hop 


M 


number of hops (number of relay nodes is M — 1) 


T m (M) 


transmissions used on hop m to achieve a successful transmission 


T(M) 


total number of transmissions required per packet, T(M) — J^T m (M) 


A 


allowance on total number of transmissions per packet, T(M) < A 


-Pout 


probability of end-to-end outage, P out = P(T(M) > A) 


C(A),C uh (A) 


random access transport capacity, and upper bound 
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M 



Fig. 2. The "capacity" and its upper bound versus the number of hops M for each M = 1, . . . , A when A = 6. Note the 
maximizing M* for both curves are close to one another. 
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Fig. 4. The optimal number of hops M*(A) that maximizes capacity C(A) and the upper bound C ub (A) versus the allowed 
number of transmission attempts A. 
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e— a = 3, SNR = 3 
O a = 3, SNR = 0.1 
=w a = 4, SNR = 3 
A ■ a = 4, SNR = 0.1 




Contention Density X 

Fig. 5. The optimal number of hops A/* normalized by \f\ vs. contention density A. It can be observed that M* = 0(VA) 
as expected, with faster convergence at high SNR. 




10 15 
Number of Hops, M 



Fig. 6. The random access transport capacity upper bound is plotted vs. the number of hops M . There is a clear optimum 
near the derived M* . 



25 




Fig. 7. Normalized random access transport capacity C ub /\/A vs. contention density A when using M*(X) hops. The V\ 
scaling law kicks in more quickly at high SNR. 



